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pq Abstract 

O Chemical exchange observed by NMR saturation transfer (CEST) or spin-lock (SL) experiments provide a MR imaging contrast 

^2k by indirect detection of exchanging protons. Determination of relative concentrations and exchange rates are commonly achieved 

(-^ |by numerical integration of the Bloch-McConnell equations. We derive an analytical solution of the Bloch-McConnell equations 

(^ that describes the magnetization of coupled spin populations under radio frequency irradiation. As CEST and off-resonant SL 

C^ [are equivalent, their steady-state magnetization and the dynamics can be predicted by the same single eigenvalue which is the 

^_^ longitudinal relaxation rate in the rotating frame Rip. For the case of slowly exchanging systems, e.g. amide protons, the saturation 

(~| pf the small proton pool is affected by transversal relaxation (R2h)- It comes out, that R2b is also significant for intermediate 

Cl^exchange, such as amine- or hydroxyl-exchange, if pools are only partially saturated. We propose a solution for Rip that includes 

n— « Rih of the exchanging pool by extending existing approaches and verify it by numerical simulations. With the appropriate projection 

^ factors we obtain an analytical solution for CEST and SL for non-zero R2 of the exchanging pool, exchange rates in the range of 1 

^ 'to 10"* Hz, Bi from 0.1 to 10 //T, arbitrary chemical-shift differences between the exchanging pools, while considering the dilution 

• by direct water saturation across the entire Z-spectra. This allows optimization of irradiation parameters and quantification of pH- 

fj dependent exchange rates and metabolite concentrations. Additionally, we propose evaluation methods that correct for concomitant 

• ^ direct saturation effects. It is shown that existing theoretical treatments for CEST are special cases of this approach. 

'Zi^^'Key words: spin-lock, magnetization transfer, Bloch-McConnell equations, chemical exchange saturation transfer, PARACEST, 
P HyperCEST 
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1. Introduction 

The relaxation of an abundant spin population is affected 
by a rare spin population owing to inter- and intramolecular 
magnetization transfer processes mediated by scalar or dipolar 
couplings or chemical exchange ll]]. As a consequence, by se- 
lective radio frequency (rf) irradiation of a coupled rare popula- 
tion not only the relaxation dynamics, but also the steady-state 
magnetization of the abundant population can be manipulated. 
Due to this preparation, the NMR signal of the abundant pop- 
ulation contains additional information on the rare population 
and its interactions. In this context, we analyze two experi- 
ments : chemical exchange saturation transfer (CEST) H and 
off-resonant spin-lock (SL). 

CEST and SL experiments are commonly applied to en- 
hance the NMR sensitivity of protons in diluted metabolites 
in vivo isl 4, 5, 6T yielding an imaging contrast for different 
pathologies |7, 8, 9, 10, LLJ. The normalized z-magnetization 
after irradiation at different frequencies, the so-called Z-spectrum, 
is affected by relaxation and irradiation parameters. In the fol- 
lowing, the large pool of water protons is called pool a and 
the pool of dilute protons pool b. To obtain a pure contrast 
that depends only on the exchanging pool b, concomitant ef- 
fects like direct water saturation or partial labeling of the ex- 



changing proton pool must be taken into account in modeling 
of Z-spectra. Similarities between CEST and SL have been no- 
ticed before Il2lll3ll . Here we consider the projection factors 
which are required for application of static and dynamic solu- 
tions derived for SL to CEST experiments and vice versa. We 
demonstrate how the experimental data have to be normalized 
that the dynamics of CEST and SL can be described by one 
single eigenvalue, namely Rip, the longitudinal relaxation rate 
in the rotating frame. A first approximation for Rip including 
chemical exchange was published by Trott and Palmer lll4ll . In 
the present article, this approach is extended by inclusion of 
R2b, the transverse relaxation rate of pool b. 

An interesting CEST effect is amide proton transfer (APT) 
of '//in the backbone of proteins, because quantitative deter- 
mination of the exchange rate may allow noninvasive pH map- 
ping ilSll . The exchange rate kh for APT is relatively small (kt, 
=28.6 + 7.4 Hz yj) compared to the transversal relaxation rate 
of the amide proton pool R2h - l/T2b- Sun et al. measured T2b 
of 8.5 ms ( R2b - 90.9 Hz) for amine protons of aqueous cre- 
atine at Bo =9.4 T. For amino protons in ammonium chloride 
dissolved in agar gel, T2b - 40 ms {R2b = 25 Hz) was found at 
Bq -3 T III6II . Thus, R2b in tissue may be in the range of or even 
surpass kb and must be taken into account for quantification of 
kb- For systems with strong hierarchy in the eigenvalues - as 
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it is the case for diluted spin populations - we present an ap- 
proximation for Rip that includes R2b and provide an analytical 
solution for CEST and SL experiments valid for exchange rates 
in the range of R2b- 

2. Theory 

CEST and SL experiments for coupled spin systems can be 
described by classical magnetization vectors M in Euclidean 
space governed by the Bloch-McConnell (BM) equations JlTll . 
We consider a system of two spin populations: pool a (abun- 
dant pool) and pool b (rare pool) in a static magnetic field Bq = 
(0, 0, Bo), with forward rate ki, and thermal equilibrium mag- 
netizations Mo^a and Mo,b, respectively. The relative popula- 
tion fraction jf^ = fh is conserved by the back exchange rate 

The 2-pool BM equations are six coupled first-order linear 
diff'erential equations 



M = A ■ M + C, 
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(?=( 0, 0, RuMo,a, 0, 0, RibMo,bf , (3) 

given in the rotating frame (x, y, z) defined by rf irradiation with 
frequency oj,f. Aw - AcOa - a>rf - cJa is the frequency offset 
relative to the Larmor frequency cja of pool a (for '// cOqIBq - 
y - 267.5 ^). The off'set of pool b Aw^ - w,/ - lOb - Aa> - 
6bCL>a is shifted by 6b (chemical shift) relative to the abundant- 
spin resonance. In contrast to Ref. |14], we allow different 
relaxation rates 7? I and 7?2 for the pools. The assumption of their 
equality is only valid if \Ria -Rib\ '^ kb or |^2n -^2il ^ kb MM- 
Longitudinal relaxation rates R\^aib - ^/Ti,a/b are in the order 
of Hz, while transverse relaxation rates R2,a/b - i/T2,a/b are 
10-100 Hz. For semisolids R2b can take values up to 10^ Hz. 
The rf irradiation field Bi = (Bi,0, 0) in the rotating frame, 
with B[ X! jjT, induces a precession of the magnetization with 
frequency a)[ = j ■ B[ around the x-axis in the order of several 
100 Hz. The population fraction fi, is assumed to be < 1 %, 
hence ka is 0.01 to 10 Hz. 

2.1. Solution of the Bloch-McConnell equations for asymmet- 
ric populations 
The BM equations ([TJ are solved in the eigenspace of the 
matrix A leading to the general solution for the magnetization 



M{t) = Yj ' 



+ M'' 



(4) 



where A„ is the nth eigenvalue with the corresponding eigen- 
vector v„ of matrix A and M" is the stationary solution. Two 
eigenvalues are real and four are complex lll4ll . They describe 



precession and, since all real parts of the eigenvalues are nega- 
tive, the decay of the magnetization towards the stationary state 
in each pool. As shown before 11911 . if Aco or wi are large com- 
pared to the relaxation rates Ri and R2 and exchange rate kb, the 
eigensystem of pool a is mainly unaff'ected. One eigenvector 
vi is closely aligned with the effective field coeff - (wi,0, Aw) 
which defines the longitudinal direction (Zeff) in the eff'ective 
frame {Xejf,yeff,Zeff) and is tilted around the y-axis by the angle 
6 - tan"'(|^) off the z-axis (Fig. [1^). Mathematical deriva- 



tion (Appendix A 1 as well as numerical evaluations (FigUJi-d) 
demonstrate that i^i and w^,^ are collinear in good approxima- 
tion if (/?2fl - /?ia) is much smaller than w^^ . 




Figure 1: (a) Geometiy of the vectors in tlie rotating frame, (b-d) Cosine of 
tlie angle jS between the eigenvector of the smallest eigenvalue and aicff. (b) 
In the far off-resonant case both vectors are parallel. Near resonance oj[ has 
to be strong to keep them parallel, (c) The assumption of coUineaiity is still 
valid if pool b with relative concentration /j < 10% is coupled to the water 
pool, (d) Large differences in R2a and i?io lead to an increasing angle between 
the vectors, but even for R2a ~ 5QHz and Si„ ~ IHz both vectors are still 
collinear in good approximation. The eigenvector and the effective field vector 
are colHnear if (iij,^ is large compared to {R2a - Ria) jAppendix A) and fi, < 
10% - both is fulfilled for CEST experiments since metabolite concentrations 
are small and frequency offsets of interest are mostly larger than several 100 
rad/s. 

The collinearity of the corresponding eigenvector and the 
effective field is the principal reason why off-resonant SL and 
CEST exhibit the same dynamics. For an appropriate analysis 
of a saturation experiment it is mandatory to identify the initial 
projections on the eigenvectors and the measured components. 
Bo and Mq are parallel to the z-axis, the preparation is a projec- 
tion of the longitudinal magnetization along z onto the effective 
frame 

M^^it = 0) = cos • M,(t = 0) = P-^j,. ■ Mo, (5) 

M,Jt = 0) = sin 6» • Mo; Myjt = 0) = 0. (6) 

The transversal components induce an oscillation decaying with 
T2p 120(1 which can be neglected in the case of small 6, by av- 
eraging over a complete cycle of cjeff, or by measuring after a 
delay of 5 • T2p. This simplification leads to the relation for the 
back projection, via P-, from Zejf to z 



M,{t) ^ cos e-M,Jt) = P,-M,Jt). 



(7) 



Since we identified the effective frame as the eigenspace of the 
magnetization, Eq. (|4| can be written as an exponential decay 
law with the eigenvalue Ai associated with the Zeff direction. 
Let the normalized magnetization be Z = -rf^ and, for the sta- 

^■" , . 

tionary solution, Z" = jf^. Then Eq. ^, taken for the Zeff 
direction, yields the dynamic solution for the z-magnetization 



Z(Aw, wi , f) = (PzPz,^. - Z") ■ e'^'-' + Z" 



(8) 



Without preparation pulses P^ - P-^^ w cos (CEST exper- 
iment). If a preparation pulse with flip angle is applied before 
and after cw irradiation the projection factors are P, - P,^^ ^ 1 
(SL experiment), hence oscillations are suppressed (Fig. |2]i, but 
still persist since Zeg is not perfectly coUinear with the eigen- 
vector. Transformation of Eq. dl} into the effective frame and 
setting jjM - yields the steady-state solution ([Appendix A|l 



Z"(Aw,wi) 



P^ • R\a -COSO 



(9) 



It is important to note that in the case where the steady-state 
is non-zero, it is locked along the corresponding eigenvector 
Equations ([8]l and (|9]l agree with the full solution previously 
found for SL by Jin et al. |5] but extend it for CEST. 

To obtain a pure dynamic quantity independent of the steady- 
state we rearrange Eq.® and define 



Z(Aw, a)i,t) = 



Z-Z' 



P,P. . - z^- 



(10) 



Eqs. ^ and (fTOl l are the central formulas in this article. 

In fact, the description of SL and CEST experiments dif- 
fers in the projection factors P, and P.^^.. The intuitive solution 
ZcEST = cos 6 ■ ZsL is valid for the steady-state, but not for the 
transient-state. If the initial magnetization M, is not fully re- 
laxed and flipped before the saturation pulse by an angle /3, P-^^ 
changes to cos(0 - P) ■ Mj/Mq. 

After understanding of the transition between the two ex- 
periments we will now solve the dynamics of CEST and SL ex- 
periments by finding the corresponding eigenvalue and verify it 
numerically. 
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Figure 2: The full numerical Bloch-McConnell solution (dots) with the pro- 
posed normalization (Eq. jlOl ) demonstrates the equivalence of two exper- 
iments; chemical exchange saturation transfer (CEST) without preparation 
pulses (a); spin-lock (SL) with preparation and measurement in the effective 
frame (b). Z of CEST undergoes oscillations because of residual transver- 
sal magnetization in the effective frame. SL shows no oscillations since the 
transversal magnetization in the effective frame is zero [AZ = 0, see text] 
(Eq. (24)). Both, SL and CEST, show the same monoexponential decay of 
the z-magnetization with Ai (Eq. )14t ) (solid red). For full BM simulations 
12111 parameters were taken from the amide proton system JSJI in brain white 
matter J2ll at Bq = 3 T: If not vaiied, R20 = 14.5 Hz, i?in = Rit, = 0.954 Hz, 
R2b = 66.6 Hz, fb = 1 %, h = 25 Hz, Si, = 3.5 ppm, Bi = 1 pT, t,„, = 1 s. 



As already demonstrated for the SL experiment 11141 . the 
eigenvalue, which corresponds to the eigenvector along the Zeff- 
axis, is the smallest eigenvalue in modulus of the system. As- 
suming that all eigenvalues of an arbitrary full-rank matrix A 
are much larger in modulus than the smallest eigenvalue, i.e. 
\Ai\ <s: \A2...n\, we obtain (see Appendix B 1 



Al ~ , 

Cl 



(11) 



where cq and ci are the coefficients of the constant and the linear 
term of the normalized characteristic polynomial, respectively. 
We derive the full solution for the smallest eigenvalue by em- 
ploying the solution of the unperturbed system {fi, = 0). The 
solution is A^j - -Reg with the decay rate in the effective frame 
Reff which was shown to be approximately 11911 



-Rejf - R\a cos" 6 + Rna sin 



(12) 



With this eigenvalue of the unperturbed system we can rescale 
the system by 



A' = A - 1 ■ i. 



f# 



(13) 



thus shifting the smallest eigenvalue by Reff. The smallest 
eigenvalue of A', still contains terms of R\a and R2a, but repre- 
sents the exchange-induced perturbation of Pip. 
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Figure 3: Hierarchy of numerically calculated BM eigenvalues (lines) of the 
standard system (see caption of Fig. |2)- (b) The rescaled matrix A' (eq. n3t ) 
has a much stronger hierarchy in the eigenvalues than matrix A (a). This im- 
proves the approximation of the smallest eigenvalue (x, Eq. II U ). 

The result is a strong hierarchy (Fig. |3) in the eigenvalues 
of A' if the coupling is small (fb «; 1). Now Eq. (fTTT i can 
be employed to calculate the eigenvalue A[ of the matrix A' to 
obtain the full solution: 



Ai - Aeff + A\. 



(14) 



Here A'^ - -c'^jc'^ is the ratio of the coefficients of the charac- 
teristic polynomial of the matrix A'. This analytical procedure 
gives us a very good approximation of the dynamics of the BM 
system. 

For further simplification we assume that relaxation of pool 
a is well described by Reg and the perturbation is dominated by 
the exchange and relaxation of pool b. We call the exchange- 
dependent relaxation rate /?„ - -Ay The eigenvalue /I I is asso- 
ciated with Zeff and is therefore an approximation of the relax- 
ation rate in the rotating frame R\p ^ -Ai given by Eqs. ( fTSI ) 
and(fT4li 



Rip(Aco) = Reff(Auj) + ReAAoj). 



(15) 



To derive a useful approximation of R^x, we neglect all relax- 
ation terms of pool a in matrix A'. Furthermore, we assume 
that Rih is much smaller than R2b and kh and therefore Rib can 
be neglected in A'. In contrast to Trott and Palmer [|l4||, we 
do not neglect R2b, but k„. By this means, the obtained eigen- 
value approximation by using Eq. (fTTT i is linearized in the small 
parameter fb giving 



nmax r^ 
RexiAcOb) = 7^ 



(16) 



with maximum value 
R'^^'^fbkb sin 



, (cOb - OJaf + ^{UJ] + AU?) + R2h(kb + R2b) 



4 



and full width at half maximum (FWHM) 



r = 2 



kb + Rib 

kb 



a)\ + {kb+R2b)^- 



(17) 



(18) 



For large \ojb - w,,! 

Rex ~ Jbkh 



w; + kbikb + R2b) 



(19) 



The wi -dependent factor yields the amount of labeling of pool 
b. Hence, we call this factor labeling efficiency, refering to Fi23ll : 



wr 



a — 



on] + kbikb + Rib) ' 



(20) 



For strong B\ and small R2b and kb,a\s approximately one 
and we obtain the full-saturation limit 



r: 



fbh = ka. 



(21) 



3. Results 



We obtained numerical values for the eigenvalues computed 
by means of the full numerical BM matrix solution Oil] and 
compared them to the proposed approximations via 



'^ex ~ 1^1, numerical I 



-R. 



(22) 



To verify equations ( 110191141171 1 the dynamics of the magne- 
tization vectors of the exchanging spin pools were simulated. 
The decay rate R^x is obtained from Z (Eq. ( fTOb ) and R^ff via 



R. 



log(Z) 



R. 



(23) 



The simulation parameters for the abundant pool were chosen 
according to published data for brain white matter ||22|| includ- 
ing a rare pool attributed to amide protons (J2|. 
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Figure 4: (a,c,e) R„- on-resonant on pool h from smallest eigenvalue in mod- 
ulus (Eq. (22)) , calculated numeiically (line) and analytically by the ap- 
proximations of Eq. )16t (squares) and the asymmetric population limit of 
Trott and Palmer 1 1411 (Eq. j251 . diamonds), x and + mark the relative error 
(1 - (R"™//?™")) when it is larger than 0.1 %. (a) For small R2i. both solutions 
for R„ agree with the numerical value; if R26 is larger than ki, the proposed solu- 
tion still matches the numerical value. The extension by R2/3 is important if the 
CEST pool is not fully saturated, which is the case for small B[ (c) or large ki, 
(e). But also for large Bi the solution, that includes R26 fits the numerical value 
with higher accuracy. (b,d,f) Imaginary parts of the numerical eigenvalues. B[ 
and R2b ranges where Eq. U6t (squares) shows deviations from the numerical 
solution con'elate with ranges where the imaginary part becomes small or even 
zero. In this case, the assumption of a strong hierarchy in the eigenvalues is not 
valid anymore. 

The proposed approximation of R^x by Eq.(fT6]l was com- 
pared to the asymmetric population solution of Ref.114] (Fig. 
|4|i. If R2b is non-zero, Re^ proposed by Eq.(fT6]l matches the 
numerical value better than the Rg^ given in Ref . 11411 (see Eq. 
dZSl) below). Especially the dependence of Re^ on B[ (Fig. |4};) 
changes by taking R2b into account. 
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Figure 5: Numerical simulations (dots) of an CEST experiment evaluated for 
dynamic Z-spectra (a, b) and steady-state Z-spectra (c, d) are in agreement with 
Eqs. UP) and (9) (solid red), respectively. Plots demonstrate high correlation 
for different Bi(0.1 - 20 jiT; a, c) and k},{\ - 2000 Hz; h, d). Deviations near 
resonance of pool a for large B\ (a) are caused by oscillations of the magneti- 
zation in the transverse plane of the effective frame. In the SL experiment these 
oscillations are suppressed (Fig.|2)- 

For an CEST experiment, the normalized numerical solu- 
tion agrees with the theory of dynamic Z-spectra (Fig.|5^, b) and 
steady-state Z-spectra (Fig. |5};, d) for different values of B \ and 
kb- The competing direct and exchange-dependent saturation 
- a central problem in proton CEST lllSi |23[ 12411 - is modeled 
correctly. Deviations in Fig. |5^ for strong B\ and Aw — > re- 
sult from transversal magnetization in the effective frame which 
was neglected before. By projection on the transverse plane of 
the effective frame using Eqs. I© we obtained the resulting de- 
viation of Z 



AZ(Aw,wi,f) - 



P-^P-^., 



P^P^.,-2^- 



■ Re(e'^-''), 



(24) 



with projections P^^^ and P, into the transverse plane of the 
effective frame and back. For MT P, — P , , - sin0. Real 
and imaginary parts of the complex eigenvalue A2 are given by 
-Rip ~ -hiRia + R\a sin^ 6 + R2a cos^ 0) 1I20I1 and Ueff, respec- 
tively. The implicit neglect of Z^ in Eq. (fTOl i is justified if 
tsat s> T2p or Pj and P^,,, are small. This can be realized ei- 
ther by SL preparation or by wi « Aw. The on-resonant case 
of CEST (0 = 90°) is not defined, because Z" in Eq. ^ and 
thus the denominators in Eqs. ( fTOt and ( l24l i vanish. Then the 
z-axis lies in the transverse plane of the effective frame and Z is 
described by Mo.a • Pe(e*"*''-''^'"'*'*). Therefore, near resonance 
SL is preferable to CEST; it also yields in general a higher SNR 
(given by the projection factors P., Pz^)- Regarding the exper- 
imental realization, CEST is simpler than SL, because Aw and 
u)\ and thus 9 can be corrected effectively after the measurement 
by Bo and B\ field mapping 1231 12511 . In contrast, SL requires 
knowledge of B\ and By during the scan for proper preparation 
or techniques that are insensitive to field inhomogeneities such 
as adiabatic pulses 11261 12711 . 
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Figure 6: Numerical i?"™ (dots), employing Eq. )23t . fit to -X'^ of Eq. )14t 
(solid blue) and to i?""' of Eq. 117) (dashed green) as a function of /j, I;,, B|, 
and S26- X and + mark the relative en'or (1 - (i?°J.'°/R™"")) when it is larger 
than 1 %. (a) As expected, the eigenvalue approximation is insufficient for 
fh > 5%. (b) Contrary to the approximation of Eg. jlll (solid red) the full 
solution follows the decrease of R"'"" with large ki, . The decrease of deviations 
for small B| (c) and high i?26 (d) may be caused by overdamping in pool b, i.e., 
eigenvalues become real which reduces the required hierarchy in the set of the 
eigenvalues. However, the deviation in S„ is too small (Fig. |4) to explain this 
deviation leading to the conclusion that other eigenvectors are contributing to 
the relaxation, (d) Inclusion of Rii, is relevant for Ry, > 100 Hz. 

The values of the rate /?™°* obtained by simulations fit well 
to the full (Eq. (fT4l i) and approximate (Eq. (fTTl i) solution for 
the observed parameters (Fig. |6). Deviations of simulation and 
analytical solution were smaller than 1% for rates varied in the 
ranges: /?ife = 0.1 - lOHz, T^i,, = 0.1 - lOHz and/?2« = 2-100 
Hz (data not shown). 



4. Discussion 

4.1. General solution 

We showed that our formalism, established by Eqs. (fT0]i.(l9]l 
and (O together with the eigenvalue approximation of Eq. (O, 
is a general solution for CEST experiments. This now allows 
us to discuss from a general point of view the techniques and 
theories proposed in the field of chemical exchange saturation 
transfer. For the SL solution this was already accomplished by 
Jinetal. iHi]. 

The proposed eigenvalue approximation assumes the case 
of asymmetric populations. This restricts its application to sys- 
tems where the water proton pool is much larger than the ex- 
changing pools - which is the case for CEST experiments. There 
are many analytical approaches for the smallest eigenvalue {R\p) 
of the BM matrix besides our approach. They use pertubation 
theory 11911 . the stochastic Liouville equation 112811 . an average 
magnetization approach |29ll. and the polynomial root finding 
algorithm of Laguerre ill 811 . The latter is even valid in the case 
of symmetric populations. However, all these treatments ne- 
glect the transverse relaxation of the exchanging pool. Since 
in CEST experiments the exchange rates are often quite small 



(e.g., kb ~ 28//z for APT), R2h cannot be neglected against kb. 
We chose therefore a simple approach which is suitable for the 
condition of asymmetric populations and took R2b into account. 
Our approach to find the eigenvalue including R2b is similar to 
that of Trott and Palmer lll4ll . However, different R\ and R2 
were allowed for the involved pools. In addition, an alterna- 
tive justification of the relation A\ - —— was obtained, which 
uses the intrinsic hierarchy of the eigenvalues ( [Appendix B[ ) 
instead of linearization of the characteristic polynomial. By 
this means, it turned out that a strong hierarchy of the eigenval- 
ues is necessary for the approximation. The hierarchy was in- 
creased by rescaling the system by the unperturbed eigenvalue 
Reff (Figl3]l. Thus the accuracy of the approximation was im- 
proved. As the parameter R2b was included and equations were 
linearized directly in the small parameter fc,,, a formula was ob- 
tained (Eq. ([15])) that differs from the asymmetric population 
limit of Ref. 11411 reading 



R\p- Reff + sin~0 



(OJb - OJa) 



2 hh 
k„+ki, 



A0.2. 



• (ka + hf 



(25) 



Equality is reached if /?2A is neglected in our approximation and 
if Eq. ( [25T l is linearized in kg. With our extension simulated 
CEST Z-spectra could be predicted well in a broad range of 
parameters. Moreover, it turned out that R2b is important if it 
is in the range of kb (Fig. [6]i). Inclusion of R2b also allows 
to model macro molecular magnetization transfer effects with 
large R2b values (Fig. [8};). 



Our solution agrees for SL with the existing treatment 1[12| . 
but only with the correct projection factors SL and CEST can be 
described by the same theory. This is contrary to the conclusion 
of Jin et al. 11211 that SL theory can be used directly to describe 
CEST experiments. The deviation is not large for small 9, but 
for a»i ai Aw the projection factors are crucial as shown in Fig. 

El 

With the correct projections the transition to CEST is straight- 
forward and provides a much broader range of validity than pre- 
vious models developed for CEST which are either appropriate 



or only for the case of 
151 1. The proposed the- 



only for small B\ (24'] or large fii ]3C 

on-resonant irradiation of pool b 123| 

ory (Eq. ([8]l) gives a model for full Z-spectra for transient and 

steady-state CEST experiments which enables analytical rather 

than numerical fitting of experimental data. 







<) Rgjj from SL-Z with correct proj. factors 
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Figure 7: The same plot of R"f-^ (Eq. (23)) as in Fig. (S};, but now for CEST 
(dots) and SL (diamonds) employing the coiTesponding projection factors in 
Eq. (To) (P- = Pj^^ = 1 for SL and P. = P^^^ = cos for CEST). Additionally, 
the result of an evaluation is shown employing the projection factors of SL for 
a CEST experiment (circles) and employing the projection factors of CEST 
for a SL experiment (squares). Only with the correct projection factors both 
experiments are described by the same theory and yield R„ (solid blue). 



4.2. Extension to other systems 



As verified for SL 111911 . the theory can be extended to «-site 
exchanging systems. By simply superimposing the exchange- 
dependent relaxation rates of several pools one obtains the Z- 
spectra for a multi-pool system. We applied this to the contrast 
agent iopamidol in water, which has two exchanging amide pro- 
ton groups 113 ill , considering a three pool system: water, amide 
proton B at 4.2 ppm and amide proton C at 5.5 ppm. Assum- 
ing for the exchange rates kc - 6 ■ kb, the superposition of 
Reff and the two corresponding /?„ yields the Z-spectrum of 
the iopamidol system (Fig. [8^). A three-pool system relevant 
for in vivo CEST studies includes water protons, amide protons 
and a macromolecular proton pool. Modeling the macromolec- 
ular pool by R'^,(R2„ = 5000 Hz, k,,, = 40 Hz) with an offset of 
-2.6 ppm and again superimposing it with R""'"''' we are able to 
model analytically Z-spectra of APT with an underlying sym- 
metric and asymmetric MT effect up to 5% relative concentra- 
tion /„, (Fig. Wp)- Hence, the model is able to describe the in 
vivo situation of several CEST pools and underlying MT com- 
peting with direct water saturation. Using the superimposed 
Rex including TJJJJ. and R™"^^ and fitting the obtained Z-spectra 
j^aimde ^^^ ^^ isolated. For macromolecular MT the extension of 
Rex by R2b is crucial, since /?2i can be as large as a; 10^ Hz. The 
implicitly assumed Lorentzian lineshape of the macromolecu- 
lar pool is only valid around the water proton resonance, for 
large offsets a super-Lorentzian lineshape must be included in 
RTxM. 
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Figure 8: Three applications of the proposed theory:(a) The system of iopami- 
dol with corresponding MTRasym evaluation, (b) The system of APT with an 
asymmetric macromolecular MTpool (for concentration fractions /„ =1%, 3%, 
5%). (c) The system of exchanging hyperpolarized xenon soluted or encapsu- 
lated in cryptophane cages, a biosensor method called HyperCEST. 



Hyperpolarized xenon spin ensembles exchanging between 
the dissolved phase and cryptophane cages (HyperCEST ex- 
periment, 1321]) can also be described by Eq. (fTOl i. Since the 
initial hyperpolarized magnetization M, is in the order up to 
10^...10*'Mo, the steady-state can be neglected for depolariza- 
tion. This yields Z ^ Z - Mj ■ g^^ip'"' in agreement with the 
result in Ref. ll33ll . Figure |8}; shows the simulated Z-spectrum 
around the cage peak in the HyperCEST experiment of a Xe- 
cryptophane system for different kt,. 

Pulsed irradiation, employed for saturation in SAR limited 
clinical scanners ll34ll . was shown to have similar effects on 
MTRasvni as cw irradiation with effective Bi 135113611 . The pre- 



^asym 



sented solution for CEST Z-spectra can therefore be used for 
optimization of pulsed saturation transfer experiments. 

4.3. Proton transfer ratio 

For a CEST experiment the parameters of particular inter- 
est are the exchange rate kt, of the metabolite proton pool and 
the relative concentration fi,. The former is often pH catalyzed 
and permits pH- weighted imaging; the latter allows molecular 
imaging with enhanced sensitivity. The ultimate method must 
allow - with high spectral selectivity - the generation of kb and 
fb maps separately and for different exchanging groups. Unfor- 
tunately, both parameters occur in the water pool BM equations 
as product, i.e. the back-exchange rate kg - kb ■ fb- There 
are some approaches which are able to separate kb and fb for 
specific cases like rotation transfer of amid protons ll37ll or the 
method of Dixon et al. (Hi applicable to PARACEST agents. 
CEST experiments are commonly evaluated to yield the proton 
transfer ratio PTR. PTR is an ideal parameter in the sense that it 
reflects the decrease of the water pool signal owing to exchange 
from a labeled exchanging pool only, thus neglecting any direct 
saturation. 



In the following, we assume one CEST pool resonance on 
the positive h.co axis. 

Employing Eq. ^ with the limit — > we obtain for PTR 
in steady-state: 



PTR = 1 - Z"(Aw) 



R\a + Rex(+^0J) 



(26) 



which yields the maximal value ^— ^ tM in the full-saturation 
Umit {Rex ~ ka ). Eq. (|26] | is consistent with PTR including the 
labeling efficiency a introduced in Ref. 12311 . 

4.4. Z-spectra evaluation - MTR and MTRasym 

Methods using asymmetry implicitly assume that the full 
width at half maximum of Rex (Aw^) is narrow compared to 
the chemical shift of the corresponding pool. This means that 
Rex (Awfe) can be neglected for the reference scan Z(-Aw) what 
is only true in the slow-exchange limit |2|l. This limit can be 
defined more generally by the width of /?„ i^^b) (Eq. (fTsTl): 



r = 2 



kh+Rn. 



2h 



h 



a)] + (kh + Ribf « |Awi - A(Da\ (27) 



This new limit depends on Bi which affects the ability to dis- 
tinguish different peaks in the Z-spectrum (Fig. ^). The limit 
is therefore a useful parameter for exchange-regime characteri- 
zation in saturation spectroscopy. 

For CEST the common evaluation parameters are the mag- 
netization transfer rate MTR(Aw) = 1 - Z(Aw) and the asym- 
metry of the Z-spectrum MTRasym(Aw) - Z(-AdLi) - Z(+Aw). 
MTRasym is generally employed to estimate PTR. Using Eq. ^ 
together with Eq. ( fTSl l we obtain for steady-state Z-spectrum 
asymmetry 

MTR^3^y„(A6>) = Z'"(-Aw) - Z'''(+Aw) 

_ (/?„(-HAw) - Rexi-i^Oj)) ■ R\aP- COS 6 (28) 
~ {Reff+Rex{-^(^))(Reff+Rex{+^Oj)) ' 

The comparison shows that MTRf^y™ yields PTR of Eq. ( l26b 
onlyife^O. 



Sunetal. UM found MTR" = Pr/?a(l-o-) which com- 



asym 



bines the labeling efficiency a found by the weak-saturation- 
pulse approximation and a spillover coefficient cr from the strong- 
saturation-pulse approximation. This formula is only valid on 
resonance of pool b in contrast to Eq. (|28] |. 

Another approach, applicable for small B\, eliminates the 
spillover effect by a probabilistic approach ll24ll . This Z-spectrum 
model taken from ll24ll yields 



PTR(Aw) 



Z"(-Aw) - Z'\+^u) 



Z'"(-^(JS) - z"(-hAw) + Z"(-Aw) ■ z'"{+^u>) 

Rexi+Aco) 



cos- 6 ■ Ria + Rexi+Aoj) 



(29) 



which turns out, after substitution of Z" by Eq.®, to be an 
approximation of PTR if 6 is small. The asymmetry normal- 
ized by the reference scan was proposed for spillover correction 



3911 . Applying eq. Q yields 

Z"(-Aw) - Z"(H-Aw) Rexi+Aoj) 



Z^'i-Acj) 



Reffi+Acj) + Rexi+Aoj) 



(30) 



which again approximates PTR if 6 is small. 

By use of Eqs. (|9]l and ^ we obtain for the asymmetry in 
transient-state a bi-exponential function 

MTRasyrn(Aw,f) =MTR^^,y„(Aw) 

+e-Rip(-^'^i' .(P_P^^^ - Z'\-Aoj)) (31) 

_g-«„(+A<.)r.(p^p_^^ - Z*-'(+Aw)). 

Neglecting direct saturation of pool a and assuming P~ - P~^^ - 1 
yields the mono-exponential approximation at the CEST reso- 
nance 11 H 



MTRasy,n(Awi = 0, = MTRr,,„(Awi = 0) ■ (1 - e 



-(R|o+*o)' 



asymv 



), 



(32) 

with the rate constant /?ip = R\a + kg. This is valid if is small, 
leading to Reg ~ R\a and, with the limit of Eq. dTTT i. Rex ~ ka 
(solid red, Fig.|6j)). 



The ratiometric analysis approach QUESTRA 14111 includes 
direct saturation and is independent of steady-state. It can be 
expressed by means of Eq. (fTOl i under the same assumptions P- 
R\a and7?„ » ka 



=P-^ =1 anAR, 



QUESTRA(0 



Z(+Aaj, t) 
Z(-Aw, t) 



„-kJ 



(33) 



Another method, pCEST Ill3ll . employs R\p in an inversion re- 
covery experiment. The pCEST signal obeys the negative of 
Eq. OD if the initial inversion is introduced by f j,^ = - cos 6. 
Hence, the full dynamics of the i?ip inversion recovery signal is 



pCEST(Aw, f) 



-MTRasym(Aw,f,P,,^ 



■ cos ff) (34) 



The pCEST signal can be positive in transient-state, but is neg- 
ative in steady-state. This inversion recovery approach was 
suggested first to increase SNR for MT effect by Mangia et 
al. [!42] and for SL akeady by Santyr et al. [43] and again 
by Jin and Kim 14411 . Their iSL signal is in our notation equal 
to Z (A(ji),(jj\,\) (Eq.®) with P-^ --\ and their projection fac- 
tors for CEST and SL are identical with P~ and P. ,. For Rex 



the approximation of Ref. 111411 is used, assuming R2b - Ria- 
Especially for the quantification employing different Bi their 
approach will benefit from our approximation of Rex- By irra- 
diation with Toggling Inversion Preparation (iTIP) Jin and Kim 
were able to remove Z" which allows for direct exponential fit 
of the difference signal of SL and iSL and thus promises re- 
duced scanning time li44ll . 

4.5. Separation for Rex 

The dependence of CEST and SL on exchange is mediated 
by Rex, the exchange-dependent relaxation rate in the rotating 



frame. Since the discussed evaluation algorithms for PTR de- 
pend on direct water saturation, we propose methods which use 
the underlying structure of the Z-spectrum and solve the solu- 
tions for Rg^. For the transient state QUESTRA can be extended 
by inclusion of /?„ and the projection factors (in Z, Eq. ( fTOl l) : 



QUESTRAr (0 



Z(+Aw, t) 
Z(-Aw, t) 



^-R„(&t^)t 



(35) 



which provides direct access to R,,,^. Even without creating Z 
one can measure the experimental Rip (Aw) decay rate and ob- 
tains Re^iAcj) a! Rip{+AcL>) - Rip{-A(D) by asymmetry analysis 
of the rate R[p (Aw). For the evaluation of steady-state mea- 
surements we suggest an extension of Eq. (l30t 



MTRr„(+Aw) 



Z"(-Aw) - Z"(+Aw) _ 
Z«(-Aw) ■ Z"(-i-Aw) ~ 
1 1 /?„(+Aw) 



Z"(-i-Aw) Z''(-Aw) cos9-P-Ri„ 



(36) 



which yields R^x in units of Ru, and is independent of spillover 
Rex can be calculated by determination of Ru, and the projec- 
tion factors . can be determined by Bi mapping and Ria 
can be measured, however Ri„ is not the same as the observed 
relaxation rate Rots in a inversion or saturation recovery ex- 
periment, especially if a macromolecular pool is present Ilia]. 
Since MTRs^^ and QUESTRAs,, evaluations employ directly 
Z-spectra data, they are useful saturation transfer evaluation 
methods for determination of R^^ with correction of direct sat- 
uration. However, they are still asymmetry-based and are not 
applicable to systems with pools with opposed resonance fre- 
quencies. In this case, the most reliable evaluation is fitting 
whole Z-spectra by using Eq.® including a superimposed 7?„ 
of the contributing pools. 

4.6. Determination of Rib, kb and fb 



As proposed by Jin et al. 11211 the width F (Eq. (fTSl l) of R^x 
(Aajb) can be used to obtain kb directly. But especially for small 
kb the extension by Rob is necessary. Fitting R'™'^ for different 
Bi y ields fb and kb separately similar to the QUESP method 
14011 and Dixons Omega Plots 13811 , but again the neglect of Rib 
in Eq. (fT9T l will distort the values for kb and f,. The width of 
Rgx is a linear function of Wj ; 



5. Conclusion 

We extended the analytical solution of the BM equations 
for SL by the relaxation rate Rib and identified the projection 
factors necessary for application of the theory to CEST exper- 
iments. Temporal evolution as well as steady-state magneti- 
zation of CEST and SL experiments can be described by one 
single model governed by the smallest eigenvalue in modu- 
lus of the BM equation system which is -Rip. Rip contains 
the exchange-dependent relaxation rate R^x. We extended R^x 
by the transversal relaxation Rib which allows application of 
the theory to slow exchange, where Rib is in the order of kb 
and not negligible. R^x of different pools can be superimposed 
to a multi-pool model even for a macromolecular MT pool. 
Compared to methods designed to estimate PTR, estimators of 
Rex are less dependent on water proton relaxation. Finally, we 
showed that determination of Rex as a function of aii and Aa» al- 
lows to determine concentration, exchange rate, and transverse 
relaxation of the exchanging pool. 

Appendix A. Eigenvector approximation 

We consider the Taylor expansion in f, of the eigenvector 
/li of the smallest eigenvalue in modulus Ai. The constant term 
of this expansion evaluated on the resonance of pool b yields 
for the components of this eigenvector in pool a 



fcoi^ 



Aco) 



f (R,„+AMR2a+M) \ 

i^(^r„'+^i) 





(A.l) 



With the approximation of Eq. ( fTSl l Ai = -Reff - Ria + (Ria - 
Ria) sin^ 6 - Ria — {Ria - Ria) COS- 6. The first component of 
(*) can be neglected if 



{R]a - Reff){Rla - Reff) 



iOl 



«; u\. 



This yields 



(Rla-Rlaf^ 



sin^ 6 cos^ ( 



jj\AiJ^ Aw^ 



(A.2) 



(A.3) 



^(-0 - 



kb + ^2, 



2h 



kb 



(cj-i) + (kb + Ribf 



and 1/7?™^ is a linear function of Wj ^ 
1 kb + Rib 



R'T(cof) fb 



(«r ) + 



fhkh 



(37) 



(38) 



Hence, also the fit of Z-spectra for different Bi yields fb, kb and 
Rib, separately. 



Since -^ > ar„ this can be reduced to the condition 



{Ria -Rial **; Wq 



(A.4) 



The second component of (*) vanishes under the same condi- 
tion ( IA.4I ). After neglect of (*) and normalization the eigen- 
vector of the smallest eigenvalue (Eq. (lA.lb ) simplifies to (Fig. 

m 



Vl 



'sine 


^cos 6^ 



(A.5) 



Along this eigenvector the Bloch-McConnell equations are one- 
dimensional 



M, ,, = ^1 ■ M, „ + C 



'<■// ^^'ff 



(A.6) 



where the constant part C-,^, is the projection of C (Eq. (O) on 
the eigenvector vi dA.SI l giving 



c,^ff = cose ■/?!„. 



(A.7) 



The solution of Eq. ( IA.6b is the combination of the general so- 
lution of the homogeneous equation (which is an exponential 
function with rate Ai) superimposed with a special solution of 
the inhomogeneous equation. The steady-state is a special so- 
lution and is obtained by setting M-,^^ = which gives 



-/ll ■ Mj;^^.^. = co&Q ■ Rxa- 



(A.8) 



By backprojection on the z-axis and normalization by Mq one 
obtains the steady-state solution Eq. (|9]i: 



Z"" = 



P-RiaCos6 



(A.9) 



Appendix B. Eigenvalue approximation 

The eigenvalues Aj of a n x «-matrix are the roots of the 
normalized characteristic polynomial and are defined by 

det(A - /I ■ I) = o r + c„_i/i""' + ... + ciA + co^Q (B.l) 

where 

CO = (-1)" det(A) = (-if A, ■ ... ■ A„ (B.2) 

and 



ci = (-1) 2j — ^ — ^ "^" ■ 



1=1 



Ai 



U '"^' 



The assumption that all eigenvalues are much larger than \Ai | « 

\A2\< ... < \An\ leads to 

ci = -co ■ (^ + ^ + ... + ^)^-T (B.4) 

Ai A2 A„ Ai 

This approximation is also valid for complex eigenvalues, be- 
cause the conjugate complex is also an eigenvalue and therefore 
J- -I- jr = —4-7T- < iji- Equations (A. 2) and (A. 4) allow general 
approximation of the smallest eigenvalue in modulus by 

Co 



^1 = -- 
ci 



(B.5) 



The error is smaller than | y,'.,l,_y, \. Justified by linearization of 
the characteristic polynomial, expression (IB. 5) was also sug- 
gested in Ref. iHi. 
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